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£h ■ Abstract 

o 

O \ In this paper, equations describing one-dimensional Non-Heisenberg 

model are studied by use of generalized coherent states in real param- 
eterization and then dissipative spin wave equation for dipole and 
quadrupole branches is obtained if there is a small linear excitation 
from the ground state. Finally, it is shown that for such exchange- 
isotropy Hamiltonians, optical branch of spin wave is non-dissipative. 

S 1 Introduction 



Many condensed matter systems fully described by use of effective continuum 
field models. Topologically nontrivial field configurations have an important 
role in modeling of systems with reduced spatial dimensionality [1]. Mag- 
netic systems are usually modeled with the help of the Heisenberg exchange 
interaction [2,3]. 

However, for spin S > 1/2, the general isotropic exchange goes beyond 
the purely Heisenberg interaction bilinear in spin operators Si and includes 
higher-order terms of the type (SiSj) n with n up to 2S [4]. Due to the 
spin states, the 2S+1 complex parameters are necessary to describe each of 
them and this corresponds with the 4S+2 degrees of freedom. Two degree of 
freedom are omitted, one because of normalization condition and the other 
for arbitrary phase decrease, hence 4S parameters are required to completely 
modeled the remainder 4S degrees of freedom of spin states. [5] 

Particularly, In case S=l with the isotropic nearest neighbor exchange 
on a lattice, is derived by use of the Hamiltonian 



h = - E J&Si+i) + k{sA+i) 2 (i) 

i 

Here Sf, Sf, Sf are the spin operators acting at a site i, J and K are re- 
spectively the bilinear (Heisenberg) and biquadratic exchange integral. The 
model (1) has been discussed recently in connection with S = 1 bosonic gases 
in optical lattices [6] and in the context of the deconfined quantum critical- 
ity [7,8]. Hamiltonian (1) is a special form presented in reference [9] and 
because of importance of quadrupole excitation in ferromagnetic Martials, it 
is considered here. This article does not consider the anti-ferromagnetic and 
nematic states. 

Considering the effects of both dipole and quadrupole branches gives a 
nonlinear approximation. If higher order multipole effects considered, the 
approximation is more accurate but at the same time, deriving the equa- 
tions are too complicated. In this paper, only effect of quadrupole branch 
for Hamiltonians described by equation (1) is considered. Study of isotropic 
and anisotropic spin Hamiltonian with non-Heisenberg terms are compli- 
cated due to quadrupole excitation dynamics[4,10,ll]. Antiferromagnetic 
property of this excitation in states near the ground, proves the existence 
of it and Dzyaloshinskii calculated the effect of this excitation [12]. Also, 
numerical calculations more accurately justify laboratory results if the effect 
of quadrupole excitation in nanoparticles Feg and Mn\i is considered[13,14]. 
In addition, this method may be promising for description of the multi-spin 
configuration of the Fe in the different ligand coordinations [15]. 

To calculate effect of quadrupole excitation, at first it is necessary to 
obtain classical equivalent of Hamiltonian (1) and then to find out the solu- 
tion of spin wave, it is necessary to analyze resulted equations in case there 
is small linear excitation from the ground states. Therefore, the stages of 
process are: 

1. Obtaining coherent states for spin s=l which are coherent states of 
SU(3) group. 

2. Calculation of average values of spin operator. 

3. Classical spin Hamiltonian equation is obtained using previously cal- 
culated values. 

4. Calculating Lagrangian equation by use of Feynman path integral over 
coherent states and then computing classical equations of motion. 

5. It is necessary to substitute resulted Hamiltonian in classical equations 
of motion to obtain nonlinear equations of magnets. Solutions of these non- 
linear equations result in soliton description of magnet that is not interested 
here. 

6. Now, ground states of magnets calculated and then nonlinear equations 
are linearized above the ground states for small linear excitation. 



7. Finally, spin wave equations and dispersion equations must be calcu- 
lated. 

In this article, we write coherent states in real parameters, because each 
parameter in this representation related to one degree of freedom. In complex 
parameterization, each parameter related to two or more degrees of freedom. 
Then in physical problems, the first representation is very helpful. 

In next sections, we develop mathematical descriptions of above stages 
and analyze these descriptions. 

2 Theory and Calculations 

In quantum mechanics, coherent states are special kinds of quantum states 
that their dynamics are very close to their corresponding classical system. 
Type of coherent state is used in problem, depends on the operators sym- 
metry in Hamiltonian. Due to the operators symmetry in Hamiltonian (1), 
For more detailed description and considering all multipolar excitation of 
coherent states we used from coherent states in SU(3) group. Vacuum state 
of this group is (1,0, 0) T and coherent state for a single site in this group is 
introduced as [10]: 

IV) = D 1 (9,<j ) )e^ §z e 2l ^ xy \o} 

(2) 

where D 1 (9,<f)) is wigner function and Q xy is quadrupole moment. Two 
angles, 9 and <j) , the Euler angles, determine the direction of classical spin 
vector in spherical coordinate system. The angle 7 determines the direction 
of quadrupole moment around the spin vector and parameter g shows change 
of magnitude of spin vector and Average value of quadrupole moment . Two 
angles <fi and 7 change between to i\ and angle 9 changes between — ir to 71. 

In order to derive Lagrangian from path integral, we drive the path inte- 
gral from the following transition amplitude: 



P(^i,ti;^,t) = (^i|exp(-i(if(ti-t)))M (3) 

Using completeness relation and doing some mathematical work, the follow- 
ing equation for transition amplitude obtained: 



P(ipi,ti;il),t) = limn-H^Y^ Dnjfy) 



—i /•*! 



xexp(-=-y Lj(9,(f),g,-f)dT) 



(4) 



In the above equation L is Lagrangian and have the following form: [10] 



L = hcos2g(cos0(f) t + j t ) - H(9, 0, g, 7) (5) 

x t = ||( x — 9,(j) ) and H(9,(j),g,j) is classical Hamiltonian. By use 
of equations (4) and (5) and the action statationary principle, the classical 
equations of motion obtained. 

When Lagrangian is obtained from path integral, another two terms ap- 
pear, one is kinetic term and has "Berry phase" properties that important 
in spin tunneling phenomena and the other depends on boundary condition 
values. In this paper, we do not consider these two terms. 

Now classical equivalent of spin vector and their products must be com- 
puted so that the classical equivalent of Hamiltonian (1) is obtained. So 
consider: 



§ = y>\s\ii>) (6) 

as classical spin vector, and also consider: 



Q'^^SiSj + SjSi-^SijI) (7) 

components of quadrupole moment. Because we can write any coherent 
state as product of single site coherent states, namely: 



m = um (8) 

i 

Then Spin operators in ground state of non-single ions Hamiltonian can 
be commuted in different lattices [10]; so 

(^|M +1 |^) = ^|^|^)^|^ +1 |^) (9) 

where \i/j) = \ip) n \ip) n+1 . 

The average spin values in SU(3) group for Coherent states (2) are defined 

as [16]: 

S+ = e i4, cos{2g)sin6 
S- = e -^cos{2g)sin9 
S z = cos(2g)cos9 

(10) 



And also 
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(11) 

In above relation, S 2 is related to dipole moment and q 2 is related to 
quadrupole moment. Classical Hamiltonian can be obtained from the average 
calculation of Hamiltonian (1) over coherent states. The classical continuous 
limit of Hamiltonian in SU(3) group is: 



H = - [ —(Jcos 2 (2g) + Kcos 4 (2g) 
J an 



a 2 



HAg 2 x sin 2 {2g){J + 2Kcos 2 (2g)) + cos 2 (2g)(J + 2Kcos 2 (2g))(6 2 x + <t> 2 x sin 2 9))) 

(12) 

Where an is length of crystal sites. The above classical Hamiltonian is sub- 
stituted in equation of motion that obtained from the Lagrangian, and the 
result is classical equations of motion: 

i 
—Q t = -a 2 cos(2g)(J + K + Kcos(Ag)) 



u 
1 

1 



x<f) xx sin6 
<f>t = a 2 cos(2g)(J + K + Kcos(4g))((f) 2 x cos6 

+6 TT csc6 
9t = 



'XX K 



COO 
1 

— lt = 4cos(2g)(J + K + Kcos(4g)) 



+ (Kcos 3 (2g)(l6g 2 x - 891-5^1 
+3<f) 2 x cos{26) - 2$ xx ) + cos(2g) 
x(8g 2 x (J -K)- 2J9 2 X + 8g 2 x Kcos(4g) 

+-J<Pl(-3 + cos(29)) - J9 xx cot0) 

+Ag xx J sin{2g) + 2g xx K(sin(2g) 
+sin(6g)))al 

(13) 

In these equations, u is ha . These equations describe nonlinear dy- 
namics of non-Heisenberg ferromagnetic chain completely. If we omitted 



quadrupole excitation (g = 0) in above equations, these equations reduced 
to Landau-Lifshitz equations. Then in comparission with Landau-Lifshitz 
equations, these equations are more complete and contain more degrees of 
freedom. Note that, solutions of these equations are different forms of mag- 
netic solitons. 

In this paper, only the linearized form of equations (13) for small exci- 
tation above the ground states is considered. To this end, at first, classical 
ground states must be calculated so in above Hamiltonian only non- derivative 
part is considered: 



H = - [ —(Jcos 2 2g + Kcos%) (14) 



To find the smallest value of the H we vary it respect to all the parameter, 
the ground state is obtained at: 



9 = 0, g = \ (15) 

In this paper, only dispersion of spin wave in neighborhood of the ground 
states is studied. For this purpose, small linear excitations from the ground 
states, as shown in eq. (15), are defined: 



2g->ir + g (16) 

In this situation, the linearized classical equations of motion are: 



—0t = -a 2 (J + 2K)<f> xx 
—<JH = a 2 (J + 2K)e xx 

1 

—9t = 

UJ 

— lt = -4(J + 2K) 

u 

(17) 
Consider functions and as plane waves to obtain dispersion equation: 



= (j) e iiu)t - kx) + (j) e- i{ult - kx) 
= o e i{u}t ~ kx) + 0oe~ i(yit ~ kx) 

(18) 



Substitution of these equations in eq. (17), result in dispersion equation 
of spin wave near the ground states. 



uj\ = (J + 2K) 2 alk 4 ul 

u 2 = -4(J + 2K)lo (19) 

It is evident from equation (19) that in addition to the dispersion acoustic 
branch, there exist non-dispersion optical branches which is related to the 
dipole and quadrapole excitations. 



3 Conclusion 

In this paper, describing equations of one-dimensional isotropic non-Heisenberg 
Hamiltonians are obtained using real-parameter coherent states. It is shown 
that both dipole and quadrupole excitations have different dispersion if there 
is small linear excitation from the ground state. 

In addition, it is shown that for isotropic ferromagnets, the magnitude 
of average quadrupole moment is constant (gt = 0) and its dynamics, is 
rotational dynamics around the classical spin vector ( / ~f t 7^ 0). 
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